Magnitude and Sign Correlations in Heartbeat Fluctuations 
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We propose an approach for analyzing signals with long-range correlations by decomposing the 
signal increment series into magnitude and sign series and analyzing their scaling properties. We 
show that signals with identical long-range correlations can exhibit different time organization for 
the magnitude and sign. We find that the magnitude series relates to the nonlinear properties of the 
original time series, while the sign series relates to the linear properties. We apply our approach to 
the heartbeat interval series and find that the magnitude series is long-range correlated, while the 
sign series is anticorrelated and that both magnitude and sign series may have clinical applications. 
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PACS numbers: 87.10.+e, 87.80.+S, 87.90+y 

A broad class of physical and biological systems ex- 
hibits complex dynamics, associated with the presence of 
many components interacting over a wide range of time 
or space scales. These often-competing interactions may 
generate an output signal with fluctuations that appear 
"noisy" and "erratic" but reveal scale-invariant structure. 
One general approach to study these systems is to ana- 
lyze the ways that such fluctuations obey scaling laws 

M 

Here, we take into account that the fluctuations in 
the dynamical output of any system can be characterized 
by their magnitude (absolute value) and their direction 
(sign). These two quantities reflect the underlying in- 
teractions in a system — the resulting "force" of these 
interactions at each moment determines the magnitude 
and the direction of the fluctuations. For an important 
representative of complex dynamics — human heartbeat 
intervals — we find unexpected results for the time or- 
dering of the heartbeat interval fluctuations by study- 
ing the scaling properties of their magnitude and sign. 
We also demonstrate that fluctuations following identi- 
cal long-range correlations can exhibit very different time 
ordering for the magnitude and sign. 

We consider the time series formed by consecutive car- 
diac interbeat intervals (Fig. |l|a) and focus on the corre- 
lations in the time increments between consecutive beats. 
This time series is of general interest, in part because it 
is the output of a complex integrated control system, in- 
cluding competing stimuli from the neuroautonomic ner- 
vous system [Q . These stimuli modulate the rhythmicity 
of the heart's intrinsic pacemaker, leading to complex 
fluctuations. Previous reports indicate that these fluctu- 
ations exhibit scale-invariant properties, and are anticor- 
related over a broad range of time scales (i.e., the power 
spectrum follows a power-law where the amplitudes of 
the higher frequencies are dominant) . 

The time series of the fluctuations in heartbeat inter- 
vals can be "decomposed" into two different time series. 
We analyze separately the time series formed by the mag- 



nitude and the sign of the increments in the time inter- 
vals between successive heartbeats (Fig. Jl|b,c). We use 
2nd order detrended fluctuation analysis [p| (and not the 
conventional power spectrum) since it has the ability to 
accurately estimate correlations in the heartbeat fluctu- 
ations even when they are masked by linear trends M. 
We find for each subject in a group of 18 healthy individ- 
uals M, that the time series of the magnitudes exhibits 
correlated behavior (Fig. ^fa) (unlike the original 
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FIG. 1. (a) An example of 2000 heartbeat (RR) intervals 
of a healthy subject during daytime, (b) The magnitude series 
of a portion of the RR series (beat numbers 800-1300) shown 
in (a). Patches of more "volatile" increments with large mag- 
nitude (beat numbers 800-1000) are followed by patches of 
less volatile increments with small magnitude (beat numbers 
1000-1300), consistent with our quantitative conclusion that 
there is correlation in the magnitude time series, (c) The sign 
series (o), as well as the ARR series (•) of a portion of the RR 
series (beat numbers 1180-1230) shown in (a). The positive 
sign (+1) represents a positive increment, while the negative 
sign (—1) represents a negative increment in the RR series 
of interbeat intervals. The tendency to alternation between 
+ 1 and —1 is consistent with our quantitative conclusion that 
there is (multiscale) anticorrelation in the sign time series. 
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FIG. 2. (a) Root mean square fluctuation, F{n), for ~ 6 
hour record (~ 32, 000 data points) for the interbeat interval 
RRi series (□) of healthy subject Here, n indicates the 
time scale (in beat numbers) over which each measure is cal- 
culated. The scaling is obtained using 2nd order detrended 
fluctuation analysis, and indicates long-range anticorrelations 
in the heartbeat interval increment series ARRi Q. As ex- 
pected, the scaling properties of the heartbeat interval incre- 
ment series remain unchanged after the Fourier phase ran- 
domization (A), (b) The root mean square fluctuation of the 
integrated magnitude series (□) indicates long-range correla- 
tions in the magnitude series | ARRi | (group average exponent 
of a - 1 = 0.74 ± 0.08 where F(n)/n cx n a_1 ). After Fourier 
phase randomization of the interbeat interval increment se- 
ries we find random behavior with exponent 0.5 (A). This 
change in the scaling (after removing the nonlinear features 
in the time series) suggests that the magnitude series carries 
information about the nonlinear properties of the heartbeat 
dynamics, (c) The root mean square fluctuation of the in- 
tegrated sign series (□) indicates anticorrelated behavior in 
sign(Ai?i?i) (group average exponent of a — 1 = 0.42 ± 0.03 
where F(n)/n cx n° _1 ). The scaling properties of the sign 
series remain unchanged after the Fourier phase randomiza- 
tion (A), which suggests that the sign series relates to linear 
properties of the heartbeat interval time series. We note the 
apparent crossovers at n m 20 beats and n ~ 100 beats. A 
gradual loss of anticorrelation in the sign series is observed at 
time scales larger than n ~ 100 beats. We note, however, that 
heartbeat increments derived from the original time series are 
anticorrelated up to scales of thousands of heartbeats. 

heartbeat increment time series, which is anticorre- 
lated, Fig. ||a). The sign series, however, exhibits an- 
ticorrelated behavior (Fig. ||c) ||. Correlation in the 
magnitude series indicates that an increment with large 
magnitude is more likely to be followed by an increment 
with large magnitude. Anticorrelation in the sign series 
indicates that a positive increment is more likely to be 
followed by a negative increment. Our result for the tem- 
poral organization of heartbeat fluctuations thus suggests 
that, under healthy conditions, a large increment in the 
positive direction is more likely to be followed by a large 
increment in the negative direction. We find that this 
empirical "rule" holds over a broad range of time scales 



from several up to hundreds of beats (Fig. ||) [|T^| . 

To show that fluctuations following an identical scaling 
law can exhibit different time ordering for the magnitude 
and sign, we perform a Fourier transform on a heartbeat 
interval increment time series, preserving the amplitudes 
of the Fourier transform but randomizing the Fourier 
phases. Then we perform an inverse Fourier transform 
to create a surrogate series. This procedure eliminates 
non-linearities, preserving only the linear features (i.e. 
two-point correlations) of the original time series [ pr| . 
The new surrogate series has the same power spectrum 
as the original heartbeat interval increment time series, 
with a scaling exponent indicating long-range anticorre- 
lations in the interbeat increments (Fig. ^a). Our anal- 
ysis of the sign time series derived from this surrogate 
signal shows scaling behavior almost identical to the one 
for the sign series from the original data (Fig. ^c). How- 
ever, the magnitude time series derived from the surro- 
gate (linearized) signal exhibits uncorrelated behavior — 
a significant change from the strongly correlated behav- 
ior observed for the original magnitude series (Fig. ||b). 
Thus, the increments in the surrogate series do not fol- 
low the empirical "rule" observed for the original heart- 
beat series, although these increments follow a scaling 
law identical to the original heartbeat increment series. 
Moreover, our results raise the interesting possibility that 
the magnitude series carries information about the non- 
linear properties of the heartbeat series, while the sign 
series relates importantly to linear properties. 

Next, we investigate the relation between the scaling 
exponent of the original series and the scaling exponents 
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FIG. 3. (a) (•) An example of anticorrelated noise (the 
scaling exponent of the increment series is 0) with (□) un- 
correlated magnitude series and (A) uncorrelated sign series 
with exponent 0.5 . (Note the sign series is anticorrelated for 
n < 20 and uncorrelated for n > 100). (b) We shuffle the 
magnitude series from (a) (gray squares) and then multiply 
its elements by the elements of the sign series from (a). The 
new surrogate series (black diamonds) is uncorrelated. 
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FIG. 4. (a) The relation between the scaling exponents 
of correlated noise, the integrated magnitude series, and the 
integrated sign series for the short range regime (n < 16). We 
generate 10 series of length 32768 with different correlations 
(input exponent) and then calculate the scaling exponents of 
the original series (•), of the integrated magnitude series (□), 
and of the integrated sign series (A). In the figure we show 
the average ± 1 standard deviation. The dashed line indi- 
cates the approximate empirical relation between the different 
scaling exponents (a sign « (original + ama g nitudc)/2). This 
empirical approximation is good for the short range regime 
only (n < 16). (b) Same as (a) for the long range regime 
(n > 64). Here the approximation does not hold anymore. In 
this regime the scaling exponents of the magnitude and sign 
series are uncorrelated no matter what is the exponent of the 
original series. 

of the magnitude and the sign series. For this purpose, 
we test our approach on well-defined signals with built-in 
correlated behavior that show uncorrelated behavior for 
the magnitude and sign. First, we consider a particular 
example of correlated noise with scaling exponent equal 
to 1, for which the increment series is anticorrelated with 
scaling exponent equal to (Fig. |j|a). Surprisingly, at 
large time scales, we find that the magnitude series and 
the sign series of the increments exhibit uncorrelated be- 
havior (scaling exponent of 0.5) although the original in- 
crement series, which is the multiplication of the elements 
of these two series, is strongly anticorrelated. Moreover, 
we find that for linear colored noise with correlation ex- 
ponent less than 1.5 (i.e., with anticorrelations for the 
increment series), the magnitude and sign series of the 
increments are uncorrelated (Fig. ^>). Next, we shuf- 
fle the magnitude series by randomly exchanging pairs 
of elements. After multiplication of the elements of the 
shuffled magnitude series with the elements of the sign se- 
ries, we find that the resulting time series is uncorrelated, 
in contrast to the original increments time series which is 
strongly anticorrelated. Note that the scaling exponents 
of the magnitude and sign series remain the same as be- 
fore the shuffling (Fig. |3|b). This test indicates that the 



correlations in a time series are not related to the cor- 
relations in the magnitude and sign series, but rather to 
the particular pairing of the elements of the magnitude 
and sign series. 

At small time scales, however, we find an empirical 
approximate relation for the scaling exponents (Fig. [i]a) , 
~ 2 ("original + ^magnitude)- We observe that for the 
heartbeat series this relation is valid over a larger range 
of scales (i.e., for time scales n < 100). 
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TABLE I. Comparison of the statistics of the root mean 
square fluctuation, F(n) (calculated using the 2nd order de- 
trended fluctuation analysis method || where n is the time 
scale in beat numbers over which each measure is calculated), 
and the scaling exponents for 18 healthy subjects and 12 sub- 
jects with heart failure ^ (obtained from 6-hour records dur- 
ing the day). The scaling features of the magnitude and sign 
change significantly for the subjects with heart failure, raising 
the possibility of bedside applications, a is the best fit to the 
range 6 < n < 1024. F(n) is estimated at the crossover posi- 
tion (n = 16) (Fig. ^b) where the largest separation between 
the two groups is estimated. Since we observe two apparent 
crossovers in the scaling behavior of the sign series, we cal- 
culate the scaling exponents in three different regions : (i) 
the short range regime for time scales 6 < n < 16 with scal- 
ing exponent, ai, (ii) the intermediate regime for time scales 
16 < n < 64 with scaling exponent, Q2, (hi) and the long 
range regime for time scales 64 < n < 1024 with scaling expo- 
nent, 03. For each measure, the group average ± 1 standard 
deviation is presented. The values which show highly signif- 
icant differences (p < 0.01 by Student's f-test) between the 
healthy and heart failure groups are indicated in boldface. We 
note, surprisingly, that the short range and the intermediate 
range scaling exponents ai and a?2 of the sign series may pro- 
vide even more robust separation between healthy and heart 
failure compared to previous reports || based on the scaling 
exponents of the original heartbeat series. 



3 



Finally, we test our analysis on a group of 12 sub- 
jects with congestive heart failure B. Compared to the 
healthy subjects, the magnitude exhibits weaker corre- 
lations with a scaling exponent closer to the exponent 
of an uncorrelated series. The change in the magnitude 
exponent for the heart failure subjects is consistent with 
a previously reported loss of nonlincarity with disease 
p2| , ^3| . The sign time series of heart failure subjects 
shows scaling behavior similar to the one observed in the 
original time series, but significantly different from the 
healthy subjects (Table [j]). 

We conclude that series with identical correlation prop- 
erties can have completely different time ordering which 
can be characterized by different scaling exponents for 
the magnitude and sign series. Moreover, we show that 
the magnitude series carries information regarding the 
nonlinear properties of the original series while the sign 
series carries information regarding the linear properties 
of the original series. The significant decrease in the 
short-range scaling exponent for the sign series in heart 
failure may be related to perturbed vagal control affect- 
ing relatively high frequency fluctuations. The decrease 
of the long-range scaling exponent for the magnitude se- 
ries of the heart failure patients indicates weaker corre- 
lations and loss of nonlinearity which may be related to 
impaired feedback mechanisms of neurohormonal cardiac 
regulation. Because information obtained by decompos- 
ing the original heartbeat time series into magnitude and 
sign time series likely reflects aspects of neuroautonomic 
regulation, this type of analysis may help address the 
challenge of developing realistic models of heart rate con- 
trol in health and disease. 
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